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nian. We demonstrate that as the length of the structure increases the system undergoes a transition 
from the ballistic to the diffusive regime. This is reflected in the calculated electron density depen- 
dencies of the conductivity and the mobility. In particular, we show that the mobility of graphene 
ribbons varies as fJ*(ri) ~ n _A , with < A < 0.5. The exponent A depends on the length of the 
system with A = 0.5 corresponding to short structures in the ballistic regime, whereas the diffusive 
regime A = (when the mobility is independent on the electron density) is reached for sufficiently 
long structures. Our results can be used for the interpretation of experimental data when the value 
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diffusive). Based on our findings we discuss available experimental results. 
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I. INTRODUCTION 

The two-dimensional allotrope of carbon - graphene 
has become a subject of intensive research since its isola- 
tion in 2004fi This is because of its fundamental signifi- 
cance, its unusual electronic properties as well as its po- 
tential for numerous applications (for a review see e.g. 2 ). 

A very interesting and not fully resolved problem is 
the impact of various mechanisms such as disorder, sub- 
strate, environment, etc., on the transport properties 
of graphene. Particular attention has been devoted to 
studies of the effect of charged impurities (located in 
the substrate or on its surface) which is widely consid- 
ered to be the main mechanism limiting the mobility 
in graphene i 3 ' ' 5 ' 6 ' Indeed, recent experiments on sus- 
pended graphene sheets have demonstrated a significant 
improvement of electrical transport in suspended devices 
compared to traditional samples where the graphene is 
supported by an insulating substrate^ 

Another important mechanism that can affect the 
transport in graphene is warping, when the graphene, 
as an elastic membrane, tends to become rippled in or- 
der to minimize the elastic energy . 10 ' 11 ! 12 ' 13 ' 14 ! 15 ' 16 ' 17 The 
warped character of a graphene surface has been proved 
in diffraction experiments^ and STM measurement ;] 16 ' 17 
for both suspended samples and samples on an insulat- 
ing substrate. It has been demonstrated that due to 
the rehybridization effects and the change in the next-to- 
nearest-neighbor hopping integrals caused by curvature 



the warping generates spatially varying potential that is 
proportional to the square of the local curvature*^ 

The transport properties of the graphene can also be 
strongly affected by its interaction with the substrate and 
other materials which may exist in its environment^ 
This includes e.g. the interaction of the graphene with 
the surface polar modes of Si02 or with water molecules 
that might reside on the surface. These interactions have 
a long-range character, and because of the corrugated 
character of the graphene and/or dielectric surfaces, the 
spatial variation of these interactions would result in a 
spatially varying effective potential affecting the trans- 
port properties of the graphene sheet. 

The purpose of the present paper is twofold. Our 
first aim is to study the effect of warping on the trans- 
port properties of graphene ribbons. The warping of the 
graphene affects both the next-neighbor and the next- 
to-nearest-neighbor hopping integrals, t and t' respec- 
tively. Previous works dealt primarily with the effect 
of the modification of the next-to-nearest-neighbor hop- 
ping integrals ! 12 ' 15 This is because the next-to-nearest- 
neighbor hopping integrals i! are are much more strongly 
affected by out-of-plane deformations in comparison to 
the nearest-neighbor integrals t. On the other hand, be- 
cause the electronic and transport properties of graphene 
are primarily determined by the next-neighbor hopping, 
it is not a priori clear which effect is dominant. In 
this work we, based on the realistic model of a warped 
graphene surface and the tight-binding p-orbital Hamil- 
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tonian, numerically study the effect of modification of the 
next-neighbor hopping integrals t on the conductance of 
the graphene ribbons. We find that the modification of 
the nearest-neighbor hopping integrals due to the out- 
of-plane deformations of the graphene surface has a neg- 
ligible effect on the conductance in comparison to the 
effect of charged impurities even for moderate strength 
and concentration. 

The second and the main aim of our study is the in- 
vestigation of the transition from the ballistic to diffu- 
sive behavior of graphene ribbons with a realistic long- 
range disordered potential. One of the motivations for 
this study are recent experiments addressing the mobil- 
ity of suspended and nonsuspended graphene devices of 
submicrometer dimensions. The dimension of these de- 
vices is smaller than the phase coherence length 1$ at 
the low temperature (1$ ~ 3 — 5 /im at 0.25 K and ~ 1 
fim at 1 K) 19 i 20 and the mean free path approaches its 
ballistic valued. This indicates that these submicrometer 
devices can be in quasi-ballistic and even ballistic trans- 
port regime requiring the Landauer approach for the de- 
scription of the transport. At the same time, the electron 
transport in these devices was analyzed in terms of the 
classical mobility \x which is appropriate for a diffusive 
transport regime. In the present study we use a realistic 
model of a disordered potential and the tight-binding p- 
orbital Hamiltonian, and perform numerical calculations 
of the conductance of graphene ribbons based on the Lan- 
dauer formalism. We demonstrate that as the size of 
the system L increases the system undergoes a transition 
from the ballistic to the diffusive regime. This is reflected 
in the calculated electron density dependency of the con- 
ductivity and the mobility. In particular, we show that 
the mobility of graphene ribbons varies as ~ n~ x , 
with < A < 0.5. The exponent A depends on the size 
of the system with A = 0.5 corresponding to short struc- 
tures in the purely ballistic regime, whereas the diffusive 
regime corresponds to A = (when the mobility is in- 
dependent on the electron density) and is reached for 
sufficiently long structures. Our results can be used for 
the interpretation of experimental data when the value 
of the parameter A can distinguish the transport regime 
of the system (i.e. ballistic, quasi-ballistic or diffusive). 

It should be noted that various aspects of the 
effect of the disorder on the electron transport 
in graphene have been extensively studied in the 
pa g tl 2L2m2i2m2L2&2&^ULa&2a We stress that the fo- 
cus of our study is the understanding of the transition 
from the ballistic to the diffusive regime when the ob- 
tained electron density dependencies for the conductivity 
and the mobility can be used to extract information on 
the character of the transport regime of the system at 
hand. Note that in contrast to many previous studies fo- 
cussing on the metal-insulator transition and the strong 
localization regime, in the present paper we consider the 
case a ribbon with many propagating channels when the 
localization length exceeds the size of the system. Note 
also that because of computational limitations we keep 
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Figure 1: The schematic sketch presenting the structure under 
consideration: (a) top view of the device region (dark gray 
area) attached to the semi-infinite graphene leads (light gray 
areas), (b) free standing, rippled graphene layer, (c) graphene 
layer supported on the substrate (white area) in the presence 
of charged impurities (dots). 



in our study the width of the ribbon W constant and 
increase its length L, such that the diffusive regime is 
achieved when L/W 1. We however expect that the 
results and conclusions presented in this paper would re- 
main valid even for bulk diffusive samples with L/W ~ 1. 



Finally, it is well established that edge disorder 
strongly affects the transport properties of graphene 
& haa 2tfi&aWWZlWWLaijiaiM. However, the 
measurements of the mobility are typically done in the 
multi-terminal Hall geometry where the edges do not play 
a role. Therefore, in the present study we consider per- 
fect edges to make sure that the electron conductance is 
influenced only by the long-range potential in the bulk. 
In our calculations we use the long-range potential cor- 
responding to remote charged impurities. We however 
demonstrate that the obtained results are not particu- 
larly sensitive to the parameters of the potential. We 
therefore can expect that our findings can be applicable 
not only to the charged impurities, but to other mecha- 
nisms discussed above (e.g. interaction with the surface 
polar modes, etc.) that can also be described by a similar 
long-range potential. 



The paper is organized as follows. In Sec. II A we 
present the basics of our computational method for cal- 
culation of the conductance and the mobility of graphene 
ribbons. The models of warping and remote impurities 
are described in Sec. II B and Sec. II C. The conductiv- 
ity and the mobility of graphene ribbons in the presence 
of warping and charged impurities are presented and dis- 
cussed in Sec. III. Section IV contains the summary and 
conclusions. 



II. MODEL 



Basics 



In order to describe transport and electronic properties 
of graphene we use the standard p-orbital tight-binding 
Hamiltonian, 



(i) 



limited to the nearest-neighbor hopping. Vi denotes the 
external potential at the site i; the summation of i runs 
over the entire lattice, while j is restricted to the sites 
next to i. We related the spatial variation of the hopping 
integral Uj with bending and stretching of the graphene 
layer due to warping as described in the next section. 

The conductance G and the electron density n are 
computed with the aid of recursive Green function 
techniqu e) 36 ' 37 . We assume that the semi-infinite leads 
are perfect graphene ribbons, and the device region is a 
rectangular graphene strip and the imperfection (warping 
and long-range impurity potential) are restricted only to 
this area (see FigJTJ) . The zero-temperature conductance 
G is given by the Landauer formula 



2e 2 



(2) 



where T is the total transmission coefficient between the 
leads. Then we calculated the conductivity 



the electron density 
n(E) 
and the mobility 



dE DOS(£), 



(3) 



(4) 



(5) 



as a functions of the Fermi energy E (W and L denote 
width and length of device, respectively). The density 
of states (DOS) was computed by averaging the local 
density of states (LDOS) over the whole device area. 
The LDOS is given by the diagonal elements of the total 
Green's functional 

All the results presented here correspond to the ribbons 
of the zigzag orientation. Previous studies do not show 
a difference of the transport properties of the zigzag and 
armchair ribbons in the presence of disorder (provided 
the disorder concentration is sufficiently high) . 29 ' 32 We 
therefore expect that all the results reported here remain 
valid for the case of the armchair orientation as well. 

The effect of warping is included in our model by mod- 
ification of hopping integrals t resulting from stretching 
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Figure 2: Definition of angles (used in Eq.©) describing 
alignment of p-orbitals of neighboring atoms, (a) (j> denotes 
the angle between direction of orbital 2 and its projection 
(dashed arrow) on the plane P spanned on the direction of 
the orbital 1 and the vector linking the centers of the orbitals 
1 and 2; (b) Top view of the plane P. All vectors depicted 
here belong to the plane P and have the same meaning as 
in (a). Dotted lines are perpendicular to the vector linking 
centers of the orbitals 1 and 2. Note that in the panel (b) the 
orbital 2 is not shown. 



(contraction) and ir — a rehybridisation. The external po- 
tential of remote impurities (with inclusion the effect of 
screening by carries) is reflected in the model by changing 
of site energies. In the next subsection we present a de- 
tailed description of models for the warping and remote 
impurities. 



B. Corrugation 

The mechanical properties of graphene can be mod- 
elled by treating this system as an elastic membrane. 
One can distinguish two modes of deformation: stretch- 
ing/contraction and bending. Both of them affect 
the strength of carbon-carbon (C-C) bonds within the 
graphene sheet by changing the distance between carbon 
atoms and the alignment of their p-orbitals, respectively. 
In the tight-binding model the hopping integral corre- 
sponding the hybridised a — it bond is given by ) 38 ' 39 

t(a, 6, 0,1?) = COs($) C0s(6>) COs(</))i2p7r,2p7r(a) 

- sin(tf) sin(6% P a,2p CT (a), 
42^,2^(0) = a T e<(l + aC+|K) 2 + ^K) 3 ), 
t 2p a,2 P a(a) = a CT e aC (-l - < - ±«) 2 + ^(aC) 3 )(6) 

where a n /t ~ -4.23, a a /t ~ -4.33 and ( » 3.07A" 1 
(t = t{ao, 0,0,0) = 2.7 eV and a = 1.42 A are the 
hopping integral and the C-C bond length for flat and 
unstained graphene ribbon, respectively). <2p7r,2p7r and 
t2pa,2 P a denote the hopping integrals for pure 7r-bonds 
(in flat graphene) and cr-bonds (for collinear p-orbitals). 
The spatial orientation of the p-orbitals is described by 
the angles 0, <fi, $ as presented in FigfJl 

The analysis of diffraction patterns of corrugated 
graphene^ provides information about the range in 
which the normal to the surface varies. The measured 



4 



a) 





Figure 3: (color online) (a) A saddle-like surface and (b) a 
local minimum created by a nonuniform contraction (dark ar- 
eas) and inflation (light areas) of a flat region, (c) developable 
surface (i.e. the surface with zero Gauss curvature). Each red 
curve and the vectors showing the radius of curvature lie in 
the same normal plane. 
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range ±5° allows to estimate from Eq.JB]) an impact of 
bending on the relative change of the hopping integral, 
t(ao,6,(f>,'d) — to) /to ~ 0.4%. The information about 
the distribution of bond lengths in corrugated graphene 
is provided by the Monte Carlo simulations M- Using this 
data we can calculate the change of hopping integrals . 
The relative change (t(a, 0,0, 0) — t )/t Q for C-C bond 
variation Aa = a — do, corresponding to the half width 
at half height (HWHH) of the bond length distribution 1 * 1 -, 
is approximately 2%. This comparison shows that the ef- 
fect of strain (leading to a change of the bond length Aa) 
has a stronger impact on the modification of the hopping 
integrals in comparison to the effect of bending (related 
to the orbital alignment). 

In a corrugated free standing membrane the bending 
and the in-plain strain are related to each other i 40 i 41 i 42 
For example, in order to produce a local mini- 
mum/maximum or a saddle- like area in a flat membrane 
it is necessary to introduce a strain (see Fig[3]). Bending 
and strain of the free standing membrane can be related 
to the Gauss curvature: 

CG(x,y) = ci(x,y)c 2 (x,y), (7) 

which is the product of the minimum c\ and maximum C2 
of normal curvatures-^ (called by principal curvatures) . 
The curvatures c±, C2 are given by the inverse radiuses 
of the local curvatures, |ci| = 1/R\, \c 2 \ = I/R2, see 
FiglH The sign of ci and C2 depends on whether an in- 
tersection of normal plane with the surface is convex or 
concave. The positive value of cg at some point cor- 
responds to the presence of a local minimum/maximun 
which requires streching of this region. For a saddle-like 
area (and negative value of cg) the local contraction of 
a flat membrane is needed. The only kind of bending 
which does not produce a strain corresponds to devel- 
opable surfaces (i.e. surfaces with zero Gauss curvature) 
such as e.g. a cylindrical tube, see Figj3](a). 

Because of a small variation of the relative change of 
the bond length Aa/ao we assume a linear dependence of 
the bond elongation / contraction Aa on the Gauss curva- 
ture: 

Aa(x,y) = kc G (x,y). (8) 



Figure 4: (color online) The model of ripples on the graphene 
surface. The black lines denote characteristic length of ripples 
8nm ; (b) the magnified small area marked in (a) by a black 
frame; M and S stand for local maxima (minima) and saddle 
points, respectively; (c) the Gauss curvature corresponding 
to a part of the surface presented in (b) - red/orange (blue) 
areas mark regions with positive (negative) Gauss curvature. 



The value of the coefficient k was chosen to reproduce the 
range of Aa variation computed in 1 - 1 -. We estimated the 
HWHH of Aa from the the C-C bond length distribution 
presented in 1 - 1 *. Then, we computed the HWHH of cg 
distribution for a large (generated as described below) 
corrugated ribbon. The coefficient k was defined as a 
ratio of the HWHHs for the Aa and cg distributions. 

In our model we do not take into account an influ- 
ence of the surface corrugation on the position of carbon 
atoms. The positions of the centres of C-C bonds in flat 
lattice were projected onto a corrugated continuous sur- 
face. In these points the Gauss curvature was computed. 
The value of cg was used to calculate the change of the 
bond lengths iJH). The angles 9, cf>, i9 (see FigfSJ) were cal- 
culated from the alignment of normals representing the ir 
orbitals. The normals were computed at the ends of the 
projected bonds. The computed angles 8, 0, t9 and C-C 
bond lengths a = ao+Aa were used to calculate the next- 
neighbor hopping integrals (J6j) in the corrugated ribbon. 

In order to model the geometry of a corrugated rib- 
bon we used the description of fluctuations of elastic 
membranes presented in 1 - 1 -. The ripples on the graphene 
surface can be expanded as a Fourier series of plane 
waves characterized by the wave vectors q. In the har- 
monic approximation 1 - 1 * the in-plane- and out-of-plane 
displacement h are decoupled. In this approximation the 
mean-square amplitude of Fourier component h q is given 
by^i: 

It can be shown that the mean-square of the out-of- 
plane displacement h scales quadratically with the linear 



5 



size of the sample^: 



(10) 



According to Eq. (|T0|) the height of the ripples increases 
quadratically with the sample size. It means that a large 
sample should be crumpled ((ft 2 ) > L) which contradicts 
the experiments. The results consistent with the experi- 
ments where large sample remains approximately flat, are 
reproduced by the Monte Carlo simulationsJ>A where the 
dependence (ft 2 ) ~ q~ A (Eq. ©) remains valid for short 
wavelengths only, and saturates for the long wavelengths 
A = 27r/<7>A*«8 nm. This mechanism is responsible 
for the existence of ripples of characteristic size ~ A* , see 
Fig. a 

Using this guidance we modeled the corrugated surface 
in the following way. The surface h(r) was generated by 
a superposition of plane waves, 



h(r) = C ^ C ?i sm(qi 



Si), 



(11) 



where r — (xi,^) is the in-plane position vector. The 
directions tp+ of the wave vectors qi = gj(cos(</5,-),sin(y>j)) 
and the phases Si were chosen randomly. The length of 
wave vectors qi covers equidistantly the range 2n/L < 
qi < 27r/(3a ), where I is a leading linear size of the 
rectangular area and ao denotes the C-C bond length (we 
assume that L 3> A*). The amplitude of the mode was 

given by the harmonic approximation C q = J 2 (ft 2 ) 

for the wave length A < A*, otherwise it was kept con- 
stant and equal to C q * , where q* — 2tt/X*. We intro- 
duced the normalisation constant C to keep the av eraged 
amplitude of the out-of-plane displacement ft = y/2 (ft 2 ) 
equal to the experimental values h « 1 nm for typical 
sizes of samples. 

In Fig. [4] the graphene surface, generated by using the 
procedure described above, is shown. Long-wave ripples 
of the size ~ A* w 8 nm discussed above are clearly seen. 
A small area marked in a box in Fig. [4^ is enlarged 
in Fig. SJd. In this region two minima (blue spots), one 
maximum (yellow spot) and one saddle region (white area 
between color spots) are shown. The relation between 
the geometry of the surface and its Gauss curvature can 
be visible in Fig. [4fc. The regions of stretching (red 
and yellow spots) and contraction (blue spot) in FigHh 
correspond to the position of minima/maximum (M) and 
saddle area (S) in FigHJa, respectively. 



C. Long-range potential 

The simplest model of charged impurities corresponds 
to the short-range (5-function scattering centres, 



N imp 

Vi = u s J2 <Hn-*v), 

i' = l 



(12) 



where Us is the strength of the individual scatterer, and 
the summation over i' runs over all impurities Ni mp in 
the device region. According to the Boltzmann transport 
theory the conductivity is a — (e 2 / 'h)(2Er / 'ft) , where r 
stands for the scattering time. For the short range poten- 
tial {12]) the scattering time r ~ l/iEpriimpUg) (where 
ni m p is the impurity concentration).— In this model the 
calculated conductivity is independent on the electron 
density n, whereas the mobility /i ~ which contra- 

dicts the experimental observations.— This indicates that 
such a simplified model might not be appropriate for a 
description of a realistic system. 

A long-range character of the electrostatic interaction 
is included in the model of bare Coulomb-like scattering 
centers, 



Vi = 



1 



47re r eo 



N imp 

E 

i'=l 



(13) 



where eo and e r stand respectively for the vacuum and 
relative permittivities. For this potential the conductiv- 
ity is proportional to electron density, a ~ n/ni mp £. 
However, the application of the bare Coulomb potential 
can be justified only for low values of the electron den- 
sity when the screening effects limiting the range of the 
potential are negligible. 

The simplest screened potential is given by the 
Thomas-Fermi approximation,— 



Vi = U TF J2 



exp ( (,tf\t - r l \) 



(14) 



where the parameters Utf and £tf describe the strength 
and the range of the scattering centres for the Thomas- 
Fermi potential. The inclusion of screening allows to 
achieve both the limits of Coulomb scattering (for low 
n) and short-range scattering (for high n). 

The singularity at r% — ry in the Thomas-Fermi po- 
tential lfl4|) can cause numerical difficulties. In our calcu- 
lation we utilize a model for screened scattering centers 
of the Gaussian shape commonly used in the literature 
where the potential on the site i read a 21 i 22 i 30 , 



Vi 



Nimp 

i'=l 



exp 



(15) 



where the height of scattering centers is uniformly dis- 
tributed in the range Uy € [— 6, 5]. The strength and 
correlation between the scattering centers is described 
by the dimensionless correlator if , 22 ' 30 



(ViVi, 



K(hv F ) 2 
2tt£ 2 



exp 



2£ 2 



(16) 



(note that (Vi) = 0). The averaging is preformed for all 
possible configurations of system which differ only in the 
distribution of the position of the impurities ry and the 
strength Uy (5, £, ni mp are kept constant). 
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Figure 5: (color online) Illustration of the long-range poten- 
tial, Eq.lfT5jl. The screening range is £ = 4oo, the impurity 
concentration is rii mv = 10 12 cm -2 . 



For this impurity potential both the screening range 
£ and the impurity strength S are independent on the 
electron density and play the role of parameters. They 
are related to K by following formula^ 2 -: 

X«40.5n mip ((5A) 2 (£/ao), (17) 

where fii mp = Ni mp /N denotes the relative concentration 
expressed via the total number of impurities Ni mp and 
the atomic sites N in a sample. In our calculation we 
use rii mp = 10 12 cm~ 2 as a typical value given in the 
experiment^. We assumed that the screening range £ = 
4ao- The strength of impurities was chosen in order to get 
typical values of correlators 2 ^ (we chose K — 1, 2, 4, 8). 



III. RESULTS AND DISCUSSION 

Let us start with the comparison of the effects of 
warping and charged impurities on the conductance of 
a graphene ribbon. Figure [6] shows the calculated con- 
ductance of a representative ribbon of the dimension 
~ 31 x 32 nm. The warping modifies the conductance 
only slightly. For low energies close to the charge neu- 
trality point the conductance steps remain practically un- 
affected. For higher energies the conductance steps be- 
come somehow distorted with the conductance plateaus 
being gradually shifted down and sharp minima appear- 
ing next to the rising edges of the plateaus. At the 
same time, in the presence of charged impurities even 
of a moderate strength the conductance of the ribbon is 
distorted substantially. The conductance steps are sig- 
nificantly washed out and the overall slope of the curve 
is lowered. Note that Fig. [6] shows the conductance of 
the ribbon for two representative impurity configurations 
with n imp = 10 12 and 4 x 10 12 cm~ 2 (dashed and dot- 
ted lines respectively). For higher energies (outside the 
first conductance step) the conductance for both cases 
shows rather similar features. However, for the energies 
corresponding to the the first conductance plateau the 



Figure 6: (color online) Effect of warping and charged impu- 
rities on the conductance of graphene nanoribbon of length 
L = 30.6 nm and width W = 31.8 nm (250 x 150 sites). The 
thin line shows the conductance steps for an ideal ribbon. The 
bold line denotes the conductance for a system with warping 
only; dashed and dotted lines refer to the conductance of a 
ribbon with impurities (without warping) with the concentra- 

-in — 2 12 2 

tion ni mp = 10 cm and 4 x 10 cm respectively. The 
strength of the impurities is £ = 4a, K = 2. The inset shows 
an enhanced view of the conductance in the vicinity of the 
first conductance step. 

conductance of the sample with rii mp — 10 12 cm" 2 (in 
contrast to the case of n imp = 4 x 10 12 cm~ 2 ) remains 
practically unaffected. This is due to the fact that in the 
zigzag ribbons the transport in the vicinity of the charge 
neutrality point takes place via edge states strongly lo- 
calized near the ribbon's boundaries. For the considered 
impurity configuration for ni mp = 10 12 cm -2 no individ- 
ual impurities are located close to the boundaries, and, 
as a result, the first conductance step remains practically 
unaffected. 

Our results demonstrate that the modification of 
the nearest-neighbor hopping integrals t resulting from 
stretching (contraction) and n — a rehybridisation has a 
little effect on the conductance in comparison to the effect 
of impurities. It should be also noted that the effect of 
the modification of the next-to-nearest hopping integrals 
t' is also weak in comparison to a realistic impurity po- 
tential. For example, Kim and Castro Netoi^ estimated 
that the variation of the effective electrochemical poten- 
tial generated by a spatial variation of t' is of the order 
of 30 meV, which is an order of magnitude smaller than 
a corresponding variation of the impurity potential^^. 
Hence, in our further analysis of the conductivity and 
the mobility of graphene ribbons we will take into to ac- 
count the impact of charged impurities only. 

Let us now turn to the investigation of transition from 
the ballistic to diffusive regime of transport taking place 
as the size of the system increases. Because of compu- 
tational limitations we keep in our study the width of 
the ribbon W constant and increase its length L, such 
that the diffusive regime is achieved when L/W ^> 1. 
We however expect that the results and conclusions pre- 
sented in this paper would remain valid even for bulk 
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Figure 7: (color online) (a) The conductivity of (a) ideal 
ballistic nanoribbons and (b) nanoribbons with impurities 
as a function of the Fermi energy E. The nanoribbon 
lengths are L = 31, 61, 123, 184, 369, 738 nm corresponding to 
250,500,1000,1500,3000,6000 sites. The conductivities in (b) 
are averaged over 10 impurities configuration. The dashed 
lines show the fit a ~ E a for the energies E > 0.2t as indi- 
cated by a vertical dotted line. (Here and hereafter we choose 
E > Q.2t because for lower energies the fitted dependencies 
deviate from the power-low behavior due to sample-specific 
fluctuations). The inset shows a dependence a = a(L). The 



impurity parameters are n%. 



10"cm" 



K = 2, £ = 4a. 



(c) The conductivity as a function of the ribbon length L 
for ballistic ribbons and ribbons with impurities. The ribbon 
width is W = 53.1 nm (250 sites). 



Figure 8: (color online) (a) The DOS as a function of the 
Fermi energy of an ideal ribbon (solid line) and of a ribbon 
of the length L = 184 nm (1500 sites) for a representative 
impurity configuration, (b) The dependence of the electron 
density n on the Fermi energy for an ideal ribbon and for 
ribbons of various length with impurities. The parameters of 
the ribbons and the impurity strength are the same as in Fig. 
[3 (Note that the electron densities for the ribbons with the 
length L > 250 sites are almost indistinguishable). 



diffusive samples with L/W ~ 1. 

Figure [7] shows the conductivity a of ideal ballistic 
graphene ribbons of different lengths L exhibiting a lin- 
ear dependence on the electron energy E. This behavior 
of the conductivity reflects the corresponding linear de- 
pendence of the conductance G of zigzag ribbons as a 
function of energy (see Appendix for details) . Note that 
the conductance of an ideal ribbon in the ballistic regime 
is independent on the length of the system. Therefore, 
for fixed energy, the conductivity a = -^G increases lin- 
early with the ribbon length L as illustrated in Fig. 0(c). 
Apparently, in the ballistic regime, the conductivity does 
not exist as a local property and can not be considered 
as material parameter because it is size-dependent. 

The conductivity of the ribbons of different lengths 
as a function of the electron energy in the presence of 
impurities is shown in Fig. [7| (b). The conductivity is 
strongly reduced in comparison to the ideal conductance 
steps and is no longer a linear function of the electron 
energy. Instead, it follows a power-law dependence 



a~E a , 



(18) 
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Figure 9: (color online) The mobility as a function of the 
electron density for (a) prefect zigzag nanoribbons and for 
(b) zigzag nanoribbons ribbon in the presence of charged im- 
purities. The inset in (b) shows a dependence A = A(L). (c) 
The conductivity a function of the electron concentration for 
nanoribbons in the presence of charged impurities. The inset 
in (c) shows a dependence 7 = 'y(L). The parameters of the 
ribbons and the impurity strength are the same as in Figs. 
[U The dashed lines in (b) and (c) show the fits /1 ~ /i A 



(as indicated by a vertical dotted line). The mobilities and 
conductivities are averaged over 10 impurity configurations. 



where the exponent a approaches a — 2 for sufficiently 
long ribbons (see inset to Fig. [3(b)). Figure [7] (c) shows 
the dependence of the conductivity a on the length of 
the system in the presence of impurities. This depen- 
dence exhibits a clear saturation of the conductivity for 
sufficiently large systems, L > 6000 sites (740 nm). In 



order to understand these features, in particular, the de- 
pendence of the conductivity on the system size L, let 
us recall that the conductance of a disordered system is 
expected to obey the scaling law 



ln(l + l/g) = L/£ loc , 



(19) 



where £; oc is the localization length and g = G/Gq is the 
dimensionless conductance (Go = 2e 2 /h being the con- 
ductance unit)^. In this study we focus on the trans- 
port regime with many transmitted channels g ^> 1. It 
follows from Eq. lfl9|) that this corresponds to the case 
£; oc S> L, i.e. the localization length exceeds the size of 
the system and the conductance is inversely proportional 
to the length, G ~ Gq/L. In this transport regime re- 
ferred to as a diffusive (or Ohmic), the conductivity a 
is therefore independent on the system size and can be 
regarded as a local quantity. According to Fig. [7] (c) this 
transport regime is achieved for sufficiently long ribbons, 
L > 6000 sites (740 nm). For smaller L the system is in 
a quasi-ballistic regime when the conductivity a depends 
on the size of the system. In this case the conductance G 
is apparently more appropriate quantity to describe the 
transport properties of the system at hand. 

Let us now turn to the analysis of the mobility of 
graphene ribbons fi. For a classical (Ohmic) conduc- 
tor the mobility is the fundamental material property 
independent on the system size and the electron den- 
sity. In contrast, for ballistic ribbons the mobility is 
size-dependent and decreases with increase of the elec- 
tron density as \i ~ n~ a5 (see Appendix). Therefore one 
can expect that the mobility of ribbons varies as 



(20) 



where the exponent A ranges from in the diffusive limit 
to 0.5 in the ballistic limit. 

In order to calculate the electron mobility fx = a/en 
we, in addition to the conductivity a, have to calculate 
the electron density n in the ribbons (note that this step 
represents the most time-consuming part of our numeri- 
cal calculations because in order to calculate n we have 
to compute the DOS for all energies < E < Ep). Fig- 
ure [8] shows the DOS for a representative impurity con- 
figuration and electron densities for ribbons of different 
lengths calculated from the DOS according to Eq. |(4|). 
For ideal ribbons, the DOS follows an overall linear de- 
pendence on the energy with the singularities correspond- 
ing to openings of new propagating channels characteris- 
tic for quasi-lD systems. Because of this linear increase 
of the DOS, the electron density for an ideal ribbon fol- 
lows a quadratic dependence on the energy, n ~ E 2 (see 
Appendix). Figure [8] (a) shows that the impurities only 
smear out the singularities in the DOS of an ideal rib- 
bon, but do not reduce the DOS. Therefore, regardless 
of the ribbon lengths the average electron density in the 
ribbons with impurities is not reduced in comparison to 
the ideal ribbons, and follows the same quadratic depen- 
dence n ~ E 2 . This behavior is expected for the transport 
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Figure 10: (color online) (a) The conductivity and (b) the 
electron density versus the Fermi energy and (c) the mobil- 
ity versus the electron density for graphene nanoribbon for 
different impurity strengths K = 1, 2, 4, 8; f = 4a, n imp = 
10 12 cm~ 2 . The ribbon's dimension is L x W = 123 x 53 
nm 2 (1000 x 250 sites). The dashed lines in (a) and (c) show 
the fits according to Eqs. (fT81 and lfT9ll for respectively 
E > 0.2* and n > 1.5 x 10 13 cm -2 . (The start of respec- 
tive fitting intervals are indicated by vertical dotted lines). 



regime at hand when the localization length exceeds the 
size of the system, £i oc ^> L, such that the averaged local 
DOS is essentially independent of the size of the system. 
Note that in the opposite regime of the strong localiza- 
tion, £i oc <C L, the electron density in ribbons is strongly 
reduced due to the effect of impurities.—. 

For low energies close to the charge neutrality point 
E = the DOS of the system with impurities shows 



sharps peaks whose positions are strongly system- 
dependent, see Fig. [8] (a). This is because in this energy 
interval the system is in the strong localization regime 
with the localization length being smaller than the sys- 
tem size (for an analysis of this transport regime see Ref. 
[H^). In this energy interval the DOS is strongly system- 
dependent because it reflects a particular configuration 
of the impurity potential. This explains some differences 
in the electron densities for different ribbons in Fig. [8] 
(b) (especially for the shortest one with L — 250 sites), 
because the integration in Eq. |(4j) includes all the ener- 
gies (below Fermi energy), including those close to E = 
where the DOS can show strong system-specific fluctua- 
tions. 

Having calculated the electron density and the conduc- 
tivity, we are in a position to discuss the mobility of the 
system at hand. The mobility of the ideal ribbons of dif- 
ferent lengths L as a function of the electron density is 
shown in Fig. [9] (a). This corresponds to the ballistic 
regime with A = 0.5. For a given electron density the 
mobility of the ballistic ribbons is proportional to the 
ribbon length, /i ~ L (because /i = a/en and a ~ L, see 
Fig- 0(c)). 

Let us now investigate how the mobility of a ribbon 
evolves as we go from the quasi-ballistic to the diffusive 
regime by increasing the system size. Figure [9] (b) shows 
the mobility of the ribbon with impurities as a function of 
the electron density. As expected, this dependence satis- 
fies Eq. (|20|) . The dependence A = X(L) is shown in the 
inset to Fig. [9](b). This dependence clearly demonstrates 
that we approach the diffusive regime (with the expected 
value of A = 0) as the length of the ribbons becomes suf- 
ficiently long (L > 740 nm). This is fully consistent with 
the behavior of the conductivity discussed above which 
exhibits a transition to the diffusive transport regime as 
the length of the ribbon increases. 

The calculated dependence A = X(L) is also consistent 
with the energy dependence of the conductivity which 
approaches the quadratic behavior a ~ E 2 in the diffu- 
sive regime (see inset to Fig. [7| (b)). Because n ~ E 2 
regardless of the regime (ballistic, quasi-ballistic, or dif- 
fusive), the mobility fi — a /en becomes independent on 
the energy (and thus on the electron density with A = 0) 
only when a ~ E 2 . To illustrate this we in Fig. [9] (c) 
present the dependence a = a(n) which is plotted by 
combining previously calculated dependencies a = <j{E) 
and n — n(E). As expected from Eqs. © and (|20)l it 
follow a dependence 

cr(n) - n 7 , (21) 

where the exponent 7 ranges from 1 in the diffusive limit 
to 0.5 in the ballistic limit. 

All the results for the conductivity, electron density 
and the mobility in the graphene ribbons presented above 
correspond to the case of one representative impurity 
strength n imp = 10 12 cm~ 2 , K = 2, £ = 4a. It is impor- 
tant to stress that the scaling laws discussed above are 
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rather insensitive to a particular realization of the poten- 
tial configuration or the impurity strength provided that 
the system is in the transport regime when the localiza- 
tion length is larger than the ribbon size. This is illus- 
trated in Fig. [10] showing the conductivity, the electron 
density and the mobility in graphene ribbons for different 
impurity strength K = 1,2,4,8. As expected, the elec- 
tron density is similar for all ribbons, especially for high 
energies, see Fig. \TU\ (b). When the impurity strength 
increases the conductivity apparently decreases, see Fig. 
flOl (a) . This decrease of the conductivity leads also to the 
decrease of the mobility as shown in Fig. \W\ (c). How- 
ever, the exponents in the scaling laws (fT8| and lfT9|) are 
not particularly sensitive to the variation of the impurity 
strength K. A small difference in scaling exponents for 
different impurity strengths has a statistical origin and 
this difference diminishes as a number of impurity config- 
urations used in calculations of each curve is increased. 

Let us now use our results to discuss available experi- 
mental data. In experiments, the electron density depen- 
dence of the mobility, [i = fx(n), and the conductivity, 
(7 = a(n), are accessible^. The dependencies lf20j) . IpTj) 
can therefore be used to extract information about the 
transport regime for the system at hand. For example, 
for the mobility the exponent A = 0.5 would correspond 
to a purely ballistic transport regime, whereas A = 
would describe a purely diffusive one. An intermediate 
exponent < A < 0.5 would indicate the quasi-ballistic 
transport regime; the more close the value of A to 0, the 
more diffusive the system is. Similar arguments applies 
to the conductivity where the corresponding exponent 7 
lies between 0.5 and 1. For example, the mobility mea- 
sured by Du et aZ.— corresponds to A = 0.5 suggesting a 
purely ballistic transport regime. The electron density 
dependence of the mobility and conductivity in graphene 
ribbons was also studied by Bolotin et The mea- 
sured conductivity follows the sublinear behavior IpTj) 
with 7 < 1. They attributed the deviation from the 
linear dependence to the effect of the short-range scat- 
tering. Taking into account that the mean free path in 
their device is comparable to the device dimension, we 
can provide an alternative interpretation of their findings 
arguing that the observed in& sublinear behavior repre- 
sents a strong evidence of the quasi-ballistic transport 
regime. 



IV. CONCLUSIONS 

In the present study we perform numerical calcula- 
tions of the conductance of graphene ribbons based on 
the Landauer formalism and the tight-binding p-orbital 
Hamiltonian including the effect of warping of graphene 
and realistic long-range impurity potential. The effect of 
warping is included in our model by modification of the 
nearest neighbor hopping integrals resulting from stretch- 
ing/contraction of the surface and the ir — a rehybridisa- 
tion. We find that the modification of the next-neighbor 



hopping due to the warping of the graphene surface has 
a negligible effect on the conductance in comparison to 
the effect charged impurities even for moderate strength 
and concentration. 

The main focus of our study is a transition from the 
ballistic to the diffusive transport regime in realistic 
graphene ribbons with long-range impurities which oc- 
curs as the size of the system increases. We keep in our 
study the width of the ribbon W constant and increase 
the ribbon length L, such that the diffusive regime is 
achieved when L/W ^ 1. We however expect that the 
results and conclusions presented in this paper would re- 
main valid even for bulk diffusive samples with L/W ~ 1, 
as soon as the localization length exceeds the size of the 
system. 

We demonstrated that the conductivity of graphene 
ribbons follows a power-law dependence a ~ E a with 
1 < a ^ 2. The case a — 1 corresponds to the ballistic 
regime whereas a — 2 corresponds to the diffusive regime 
which is reached for sufficiently long ribbons. In the 
ballistic regime the conductivity scales linearly with the 
length of the system L, whereas in the diffusive regime 
the conductivity saturates with L. 

We find that the average electron density in the ribbons 
with impurities is practically not reduced in compari- 
son to the ideal ribbons, and follows the same quadratic 
dependence n ~ E 2 regardless of the transport regime 
(ballistic, quasi-ballistic or diffusive). This behavior is 
consistent with the exponent a — 2 reached in the dif- 
fusive case, because in this case the mobility /1 = a /en 
becomes independent on the energy (and hence on the 
electron density) as expected for the diffusive regime. 

In experiments the electron density dependence of the 
mobility, /i = /i(n), is accessible. We find that the mo- 
bility of graphene ribbons varies as /x(n) ~ n~ x , with 
< A < 0.5. The exponent A depends on the size of the 
system with A = 0.5 corresponding to short ribbons in 
the ballistic regime, whereas the diffusive regime A = 
(when the mobility is independent on the electron den- 
sity) is reached for sufficiently long ribbons. Our results 
can be used for the interpretation of the experimental 
data when the value of the parameter A can be used 
to distinguish the transport regime of the system (i.e. 
ballistic, quasi-ballistic or diffusive). The correspond- 
ing electron density dependence for the conductivity is, 
a{n) r~j rT 1 , where the exponent 7 ranges from 1 in the 
diffusive limit to 0.5 in the ballistic limit. Based on our 
findings we discuss the available experiments and pro- 
vide an alternative interpretation of some experimental 
conclusions.— £ 

Our calculations also demonstrate that in the quasi- 
ballistic regime (which corresponds to many experimen- 
tal studies) the mobility and the conductivity of the 
structure at hand strongly depend on the system size. 
Therefore in this regime the conductivity does not exist 
as a local property and the mobility can not be consid- 
ered as a well-defined material parameter because of its 
dependence on the system size. 
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Appendix A: ELECTRON CONDUCTIVITY, 
MOBILITY AND ELECTRON DENSITY IN THE 
BALLISTIC REGIME 



In the low-energy limit close to the charge neutrality 
point E = the electron density of the zigzag graphene 
nanoribbon of the width N reads 41 



n(E) = 



1 



VE^-El6(\E\-\E m \), (Al) 



7T\/3 ^0«0 

where the threshold energies have the form, 
3tt 1 

E m = —toao^im + ±), m = 0, 1,2, . . 
8 W 2 



(A2) 



Noticing that the function ^JE 2 — E^ approaches E for 
E ^> E m we can write the electron density in the approx- 
imate form: 



i(E) 



\E\ - \E„ 



(A3) 



Summation of ^-functions in Eq. (| A3[) gives a number 
of propagating modes at the given energy. Expressing 
this number with energy by making use of Eq. (|A2|) and 
approximating E ~ E m we obtain, 



n(E) 



.32 



1 



7T 2 3-s/3<0 O 



WE 2 



(A4) 



A comparison of the approximate expression (j A4|) with 
the exact one (|A1[) is shown in Fig. fTTIfa). 

In the low-energy limit close to the charge neutral- 
ity point E = the conductivity of the zigzag graphene 
nanoribbon of the width N reads,— 



a(E) 



~hW 



26{\E\ - \E m \) - 1 



Using similar approximations as above we obtain, 

ji- LB. 



(A5) 



(A6) 



a) 



n(E) 




b) E « E, E 2 



a(E) 
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Figure 11: (Color online) (a) The electron density, (b) the 
conductivity (c) and the mobility as a functions of the Fermi 
energy for a zigzag nanoribbon. Red lines correspond to the 
approximate expressions (( A4[) . (JA6]| , (( AT[) . whereas the black 
lines to the exact expressions l|Aip . (fA5j. 



Finally, substituting Eq. (|A6|) and (|A4|) into the defini- 
tion of the mobility we obtain, 



en 



a e r- L i e 
— w -nV3t a — E = - 



AV2 L 

7 

h ^3 VW 



-0.5 



(A7) 



A comparison of the approximate expression for the con- 
ductivity and the mobility with the exact ones is shown 
in Fig. [TI](b), (c). 

It should be noted that the expressions l|Al|l . (j A3fl . 
(| A4|) for the electron density n do not include a contribu- 
tion from the edge states existing in the zigzag nanorib- 
bons for the energies close to E = 0. Because of this the 
electron density of the nanoribbon shown in Fig. Q1] van- 
ishes at E = and the mobility /U exhibits a singularity. 
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Accounting for the contribution from the edge state in away the charge neutrality point this contribution does 
exact numerical calculations leads to the finite values of not practically affect n and /i, which justifies the utiliza- 
n and fi at E — 0. However, for high energies sufficiently tion of expressions (|Al[) - <|A7[) . 
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